M>-AM6  6*9  CLARKSON  COLL  OF  TECHNOLOGY  POTSDAM  N  T  DEPT  OF  MECH — ETC  F/6  11/9 

OPTIMALITY  CRITERION  TECHNIOUES  APPLIED  TO  STRUCTURES  COrf*OSED  — ETC«U> 
1900  M  t  SYEO»  K  D  WILLMERT *  M  R  KHAN  NOO0U-7O-C-OOO* 

UNCLASSIFIED  MlE-059  ML 


TH 

CO 

CO 

00 

o 


OPTIMALITY  CRITERION  TECHNIQUES 

APPLIED  TO  STRUCTURES 

COMPOSED  OF 
DIFFERENT  ELEMENT  TYPES 

M.  I.  Syed 
K.  D.  Willmert 
M.  R.  Khan 


Department  of 
Mechanical  and  Industrial  Engineering 

Clarkson  College 
Potsdam ,  N.Y. 

Office  of  Naval  Research 
Contract  No.  N00014-76-C-0064 


DT1C 

SELEC  I\ 

JUL7c 31 


Report  No.  MIE-059 
May  1980 


AOX 
/P?/  f 

I ••  c.. 

5(5 

°\l  "1 

V\\  \ 


r*  8121 

:>  »?,  Hi 

’  /i?/ 

M/ 


80  7  1  388 


Clarkson  College  of  Technology 
Potsdam,  New  York  13676 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Navel  Research 
Room  303  Federal  Building 

Rochester,  New  York  14614  _ 


14.  MONITORING  AGENCY  NAME  «  ADORESSfl/  dltloronl  from  Controlling  Ollleo) 

Office  of  Navel  Research 
Structural  Mechanics 
Department  of  the  Navy 

Arlington,  Virginia  22217  _ 


16.  DISTRIBUTION  STATEMENT  (ot  this  Roport) 


Mav  1980 


IS.  NUMBER  OF  RAGES 

7 


IS.  SECURITY  CLASS,  (ot  Ml*  roport) 

Unclassified 


IS*.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


This  document  has  been  approved  for  public  release  and  sale;  distribution 
unlimited . 


17.  DISTRIBUTION  STATEMENT  (ol  tho  obotroc t  ontorod  In  Block  20,  It  dltloront  from  Roport) 


ts.  supplementary  notes 


Paper  presented  at  the  AIAA/ASME  21  st  Structures,  Structural  Dynamics  and 
Materials  Conference,  Seattle,  Washington,  May  12-14,  1980.  Copied  with 
permission  of  AIAA. 


19.  KEY  WORDS  (Conti, 


mmty  and  Idantlty  by  bloc m  numbar) 


Optimization 

Optimal  Structural  Design 
Optimality  Criterion 


S,  "ABSTRACT  (Continue  on  roeoroe  old*  II  nocoooory  and  Identity  by  block  number; 

The  optimality  criterion  design  method,  developed  by  Khan,  Wlllmert  and 
Thornton  which  exploits  the  concept  of  one  most  critical  constraint,  has  been 
extended  in  this  paper  to  handle  more  complex  structures,  specifically  ones 
composed-x>f  several  different  types  of* elements.  The  elements  considered  were 
the  <£russ)  and  two  different  plane  stress* elements:  a  constant  strain 
triangular  element  and  a  rectangular  symmetric  ''-shear  panei.  The  optimality 
criterion  method  treats  multiple  load  conditions  and  stress  and  displacement  ^ 


nn  FORM 

VV  I  JAN  73 


EDITION  OF  I  NOV  63  IS  OBSOLETE 

S/N  0102-  LF-  014-  660) 


_ Unclassified  _ML 

SECURITY  CLASSIFICATION  OR  THIS  BAOt  fBSlwi 


SECURITY  CLASSIFICATION  OF  this  page  fWh  an  Data  Bnfnd) 


constraints.  A  wing-box  structure  is  considered  as  an  example  to  demonstrate 
the  efficiency  and  simplicity  of  the  method,  l 


Accession  For 

Giuil 

DL'C  ZiM 

UnaattvunMi 

Juscii'iaation 


0  *.  t  y  y t,  ^  '•*/ 

Avt'ilnbii  If  y_  Codc.s 
I  Avail  and/or 
Diat.  special 


S/N  010?-  LF-  014-  6601 


OPTIMALITY  CRITERION  TECHNIQUES 
APPLIED  TO  STRUCTURES  COMPOSED 
OP  DIFFERENT  ELEMENT  TYPES 


M.  I.  Syed* 

Moog  Incorporated 
East  Aurora,  New  York 

K.  D.  Willmert**  and  M.  R.  Khan*** 
Clarkson  College  of  Technology 
Potsdam,  New  York 


Abstract 


II .  Theory 


The  optimality  criterion  design  method,  devel¬ 
oped  by  Khan,  Willmert  and  Thornton  which  exploits 
the  concept  of  one  most  critical  constraint,  has 
been  extended  in  this  paper  to  handle  more  complex 
structures,  specifically  ones  composed  of  several 
different  types  of  elements.  The  elements  consid¬ 
ered  were  the  truss,  and  two  different  plane  stress 
elements:  a  constant  strain  triangular  element  and 
a  rectangular  symmetric  shear  panel.  The  optimality 
criterion  method  treats  multiple  load  conditions 
and  stress  and  displacement  constraints.  A  wing- 
box  structure  is  considered  as  an  example  to  demon¬ 
strate  the  efficiency  and  simplicity  of  the  method. 


I.  Introduction 


The  application  of  optimality  criterion  tech¬ 
niques  to  structural  design  has  been  of  great  inter¬ 
est  to  many  researchers  during  the  last  few 
years. (See  Venkayya9  for  an  excellent  review 
article  in  this  area.)  Khan,  willmert,  and  Thornton 
have  developed  several  very  efficient  techniques 
based  on  optimality  conditions  for  the  design  of 
both  two  and  throe  dimensional  trusses,  frames,  and 
even  high  speed  mechanical  mechanisms. 10-14  Numer¬ 
ous  constraints  have  been  considered  including 
limits  on  stress  and  displacement  (and  combinations 
of  these),  natural  frequencies,  and  upper  and  lower 
bounds  on  the  design  variables.  Recently  problems 
involving  more  than  one  design  variable  per  member 
have  also  been  considered.15  The  techniques,  based 
on  the  assumption  that  only  one  constraint  is  most 
critical  at  any  stage  of  the  process,  have  proven 
to  be  very  efficient  in  terms  of  computational  time 
and  core  storage  requirements.  It  has  been  shown 
through  numerous  examples  that  the  optimal  designs 
obtained  by  these  new  methods  wore  nearly  identical 
(sometimes  better)  to  those  obtained  by  previous 
methods. 

In  all  of  the  problems  considered  thus  far, 
however,  the  structures  were  composed  of  a  single 
clement  type,  that  is,  either  truss  or  frame  ele¬ 
ments,  or  in  the  case  of  mechanical  mechanisms,  a 
special  vibrational  element.  The  purpose  of  the 
research  presented  here  was  to  develop  an  optimal¬ 
ity  criterion  technique,  based  on  the  efficient 
method  ol  Khan,  Willmert,  and  Thornton,  capable  of 
designing  structures  com|«sed  of  several  different 
clement  types. 
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The  type  of  structures  considered  in  this  work 
are  ones  composed  of  truss  elements  and  two  differ¬ 
ent  plane  stress  elements,  particularly  a  constant 
strain  triangular  element  and  a  symmetric  rectangu¬ 
lar  shear  panel.  The  cross  sectional  areas  Ai  of 
the  truss  elements,  shown  in  Figure  1,  and  the  thick¬ 
nesses  Tj  of  the  plane  stress  elements,  shown  in 
Figures  2  and  3,  were  used  as  the  design  variables. 
The  objective  was  to  minimise  the  total  weight  of 
the  structure: 

Tn  J 

W  =  p  £  A.L  +  £  T.S 

[i*=l  j«l  3 

where  Lp  and  Sj  are  the  length  of  the  truss  element 
i  and  surface  area  of  the  plane  stress  element  j 
respectively.  N  and  J  are  the  total  number  of  truss 
and  plane  stress  elements  respectively. 

Both  stress  and  displacement  limitations  were 
imposed  on  the  design.  The  stress  limitations  on 
the  truss  elements  were: 


“ik 


i  -  1,...,N 
for  all  k 


(2) 


where  Op  is  the  allowable  stress  for  element  i, 
and  apit  is  the  actual  stress  in  element  i  under  load 
condition  k.  For  the  plane  stress  elements  the  Von 
Miscs  equivalent  stress  was  used: 


2  2  1/2  — 
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j  ■  1 . J 

for  all  k 


where  oejk  is  tlle  equivalent  stress  for  element  j 
under  load  consition  k.  The  nodal  displacements  umk 
of  the  structure  were  limited  by: 
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m  «  1, . . . ,M 
for  all  k 


(4) 


where  M  is  the  total  number  of  displacement  con¬ 
straints.  Lower  bounds  wire  also  placed  on  the  areas 
and  thicknesses. 


The  constant  strain  triangular  element  used  is 
shown  in  Figure  2.  It  is  assumed  that  this  element 
has  constant  strain  in  the  field  given  by  the  equa¬ 
tions: 
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whore  Up,  ’/p,  Uq,  . ..,  Vr  are  the  di ^placements  in 
the  local  X  and  Y  directions  of  the  nodes  I',  Q,  and 
R,  and  b,  h,  and  s  are  the  dimensions  of  the  element 
as  shown  in  Figure  2.  The  other  plane  stress  ele¬ 
ment  used  la  this  research  was  the  symmetric  rec¬ 
tangular  shear  panel  element  shown  in  Figure  3. 

This  element  is  assumed  to  be  symmetric  with  respect 
to  the  X  axis,  with  a  stress  distribution  of  the 
form: 


:‘lY  ^  °2 


Figure  2:  Constant  Strain  Triangular  Element 


where  Dj,  D2,  and  D3  are  constants.  It  is  also 
assumed  that  the  displacement  in  the  X  direction  on 
this  X  axis  is  zero,  i.e. 

u (X, 0 )  =  0  (7) 

Using  the  above  assumptions,  the  element  sti ff- 
ncss  matrices  can  tie  developed  for  all  three  ele¬ 
ments;  from  which  the  displacements  and  stresses 
can  be  determined  which  arc  required  for  the  con¬ 
straints  in  the  optimization. 


figure  1:  Truss  Element  With  Uniform  Cross- 
Sectional  Area 


Figure  3:  Symmetric  Shear  Panel  Element 


Design  Algorithm 


Khan,  Willmert  and  Thornton*0  derived  recursive 
relations  for  the  design  of  structures  composed  of 
truss  or  frame  elements  based  on  the  concept  of  one 
most  critical  constraint  at  any  iteration.  From 
these  they  developed  an  algorithm  for  the  design  of 
trusses  or  frames  under  multiple  stress  and  displace¬ 
ment  limitations.  In  their  work,  the  cross-sectiona 1 
areas  were  the  design  variables.  In  the  research 
presented  here,  the  areas  of  the  truss  elements  and 
thicknesses  of  the  plane  stress  elements  were  used. 
Thus  recursion  relations  for  the  thicknesses  were 
derived  in  this  research.  These  were  combined  with 
those  developed  by  Khan,  Willmert  and  Thornton  to 
obtain  the  following  design  algorithm. 

1.  Choose  any  design,  specifically  values  for 
all  the  design  variables  namely  areas  Ai,  i«l, 
2.....N  and  thicknesses  Tj ,  j»l,2,...,J,  such 
that  all  the  areas  are  the  same  and  similarly 
for  the  thicknesses.  Choose  a  value  of  the 
relaxation  parameter  (n) .  (Say  0.05  to  0.15). 


2 


2.  Analyse  the  design  for  each  load  condition. 


3.  Chock  displacements  in  each  load  condition 
at  those  nodes  where  displacement  limitations 
are  imposed  and  determine  the  node  and  direc¬ 
tion  for  which  the  c.ilculafeil  <lis|.laceinent 
most  close iy  approaches  (or  exceeds)  the  allow¬ 
able  displacement.  This  is  tile  most  critical 
displacement  (upq) . 

4.  Knowing  the  magnitude  of  the  most  critical 
displacement  (Up,-)  from  Step  3  and  the  value 
of  the  allowable  displacement  (Up),  scale  the 
chosen  design  for  all  the  design  variables  so 
that  the  most  critical  displacement  becomes 
active.  All  other  displacement  constraints 
will  then  be  inactive.  Let  the  scaled  design 

be  denoted  by  A!  and  T(  where 
i  3 


i  =  1,2, ...,N  (8) 


j  =  1,2 . J  (9) 


If  the  structure  was  analyzed  with  the  scaled 
design  then  the  displacement  vectors  calculat¬ 
ed  at  Step  2  would  have  been: 


A1  =  — B3-  A 
i  u  i 
P 

and 


T*  =  t 

j  u„  3 


-  UP  - 

u'  =  — £ —  u  for  all  k  (10) 

k  iu  |  k 

I  pql 


and  stiffness  matrix  from  the  scaled  design 
would  have  been: 

|u  | 

[K'l  =  -^-[Kl  (11) 

P 

-► 

5.  From  the  scaled  displacements  uj^  and  design 
variables  A^  and  T*,  compute  the  maximum  stress 
mJxlai]c|  in  each  truss  member  i,  and  rajjx  °ejk 
for  each  plane  stress  element  j.  Also  deter¬ 
mine  the  stress  response  ratio  (stress  obtained 
at  Step  2  divided  by  the  limiting  stress)  for 
each  member  and  lot  the  most  critical  response 
ratio  be  obtained  for  the  nth  member.  This  is 
denoted  by  R„. 


If  R  >1  comput  e  V  -K 
n  In 


N  .3 

j:  a’i..+  r.  -r’.s. 
i=!  1  1  j  =  l  ^ 


If  R  <1  compute  V  = 
n—  1 


N  J 

l  a!l.  *  T'.s . 
Li-1  1  1  3=1  3  u 


(12) 


6.  Using  the  sealed  design,  apply  a  unit  load 
only  at  the  node  and  in  the  direction  of  the 
active  displacement  constraint.  Let  the  set 
of  resulting  nodal  displacements  be  denoted  by 
u|.  Note  that  this  is  the  only  unit  load  that 
must  be  applied,  and  that,  the  structural  stiff¬ 
ness  matrix  inverted  at  Step  2  is  used  here  as 
scaled  in  Step  4  to  compute  uj . 

7.  Compute  the  displacement  derivatives  for 
area  design  variables  as  follows: 


&U 

—13., 

3a: 


*♦  T 

ulg  gij 

Ai 


i  •  1,2, ...,N 


(13) 


where  u(  is  a  vector  of  the  actual  displace¬ 
ments  for  member  i,  and  [K^l  is  the  scaled 
element  stiffness  matrix  for  element  i.  Com¬ 
pute  displacement  derivatives  for  thickness 
design  variables  using: 


3u 

_E3. 

st: 
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ISLS 
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1,2,... ,J 


(14) 


Also  the  Lagrange  multiplier  associated  with 
tho  critical  displacement  for  area  and  thick¬ 
ness  design  variables  is  calculated  as  follows: 
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pq 


i=l 


L.  + 
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J 
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j=i 


T*  S 
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(15) 


8.  Group  the  memebers  as  follows: 

i.  If  3upq/3A^  >  0  or  a;  >  Oj,  then  truss 
member  i  belongs  to  group  G^.  If  3Upq/3Tj  _> 
0  or  o6j  >  Oj  then  plane  stress  member  j 
belongs  to  Gj . 


ii.  Otherwise,  member  i  or  j  belongs  to 
group  G2.  Note  that  either  group  could  be 
empty  and  a  particular  member  would  belong 
to  only  one  group  at  a  time. 


9.  Use  the  stress  ratio  formula  to  resize  the 
elements  of  group  G^.  Resizing  of  the  area 
design  variables  is  accomplished  as  follows: 


^i^v+l 


max  0 


ik1 


1,2, ...,N 


(16) 


and  resizing  of  the  thickness  design  variables: 


IT  } 

1  jJv+l 


Piax  a 
k  ejk 


o . 
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j  -  1.2, 


(17) 


10.  Resize  the  elements  of  group  G2  for  the 
area  design  variables  as  follows: 


(Vvl 


(-3u  /3AjfJn 

A  - & - i-  (A! 

pq  L  i 

1  Jv 


>..  i  "  1.2, . 


..,N 

(18) 


and  resize  the  elements  of  group  G2  for  thick¬ 
ness  design  variables  as  follows: 

h 

(T’)v  j  -  1,2 . J 

v  (19) 

11.  Compute  the  volume  from  areas  and  thick- 
nosses^btained  in  Step  9  t  10.  If  the  quantity 
ItVi-v^J/vJI  is  less  than  e  (a  small  number 
ranging  between  .001  to  .050)  then  stop,  this 
design  is  the  optimum  design)  otherwise  go  to 
Step  2  and  repeat  the  process. 


(tjW 


<-3u  /3TM 
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pq 


IV.  Ilxamj'l  »•*. 

To  show  th«?  ef  fi  cr  i  v«*n<‘«is  and  ffficioncy  of  th? 
a  I <jur i  t  hm,  a  win-j  Uox  struct  ur«*  as  shown  in 
P i 1 1 1 j r * •  4  was.  run  .  1 .  This  a  rurt  «ir»*  was  assum¬ 

ed  to  Lm>  syini'K't  r with  i*  »jm  «  t  to  tin*  iniddl*-  sur- 
fvK’e.  As  a  result  only  t!u*  uj*i>*  r  halt'  of  tho  struc¬ 
ture  was  considered  in  ion.  In  this  half 

there  were  5  truss  elements,  5  constant  strain  tri¬ 
angular  elements  and  8  rectangular  shear  panel  ele¬ 
ment.  The  stress  limit  used  for  all  elements  was 
10,300  psi,  lower  bounds  on  the  truss  areas  was  0.1 
in2  and  on  the  {lane  stress  element  thicknesses  was 
O.w.*  in.  'jvo  diH.-rent  load  con  til  i«>na  wore  c;on- 
side  red.  These  were  loads  of  OooO  lbs.  at  node  7 
(as  shown  in  Figure  4),  and  10, out)  lbs.  at  node  5, 
both  in  the  Z  direction.  The  problem  was  solved 
three  different  times  wirh  three  different  di splace- 
ment  limits.  These  were  2.0,  1.0  and  4.0  inches 
at  all  nodes  in  all  directions. 

O  truss  element 

A  CST  ELEMENT 

□  SSP  ELEMENT 


p, 


Figure  4:  Finite  Element  Model  For  Eighteen 
Element  Wing  Box  Structure 

The  2.0  inch  displacement  limit  case  has  been 
solved  by  other  researchers  including  Gallatly  and 
Herke^,  Gailatly*7,  and  Schmit  and  Miura^s.  The 
final  designs  obtained  are  shown  in  Tabic  1.  It  is 
difficult  to  compare  these  results  directly  since 
different  finite  element  models  of  the  structure 
were  used  in  each  case,  and  thus  some  variation  in 
optimal  designs  arc  observed.  Although  Schmit  and 
Miura’s  model  I  was  id'-ntieal  to  that  used  in  the 
research  presented  hem.  The  optimal  structures 
are  similar  however.  The  greatest  variation  among 
the  designs  shown  appears  in  the  area  of  truss  ele¬ 
ment  1.  The  optimal  design  obtained  l>v  the  techni¬ 
que  of  this  res-  at  eh,  although  slightly  higher  in 
weight,  is  similar  to  the  other  designs. 

In  order  to  provide  a  means  of  direct ly  com¬ 
paring  results,  this  example  problem  was  also  solved 
using  two  other  methods  in  addition  to  the  techni¬ 
que  presented  in  section  III  of  this  paper.  These 
were  the  SUMT  technique  using  1'owi •  1 1  ' s  method^, 
and  by  the  stress  ratio  method  with  displacement 


scaling.  In  the  latter  technique,  the  stress  ratios 
alone,  equations  (16)  and  (17),  were  used  to  change 
the  design  variables  from  one  iteration  to  the  next, 
hut  after  each  iteration  the  design  was  scaled  so 
that  tile  displacement  constraints  were  all  satisfied. 
Table  2  presents  the  resulting  optimal  designs  for 
an  allowable  displacement  of  2.0  Inches  obtained  by 
the  three  methods  of  this  research.  Two  different 
values  of  the  relaxation  parameter  n  were  used  with 
the  algorithm  of  section  111  as  shown  in  this  table. 

The  optimal  design  obtained  fc.  TowelVs  method 
was  almost  identical  to  that  reported  by  Schmit  and 
Miura  using  their  Model  Ij  however,  the  technique 
required  4605  analyses  to  obtain  this  design.  The 
stress  ratio  method  with  displacement  scaling  also 
gave  a  similar  design  with  only  slightly  higher 
weight.  It  is  noted  that  the  final  design  obtained 
by  this  method  was  not  fully  stressed,  because  of 
the  displacement  scaling  present.  The  optimal  de¬ 
signs  generated  by  the  technique  of  section  111  of 
this  paper  for  two  different  relaxation  parameters 
had  a  slightly  higher  weight  yet.  This  seems  to 
indicate  that  the  stress  ratio  method  with  displace¬ 
ment  scaling  is  perhaps  the  best  method  to  use  in 
this  case,  since  it  is  very  simple  and  is  able  to 
locate  very  nearly  the  optimal  design  in  a  small 
number  of  iterations. 

The  wing  box  problem  was  also  solved  under  a 
reduced  allowable  displacement  of  1.0  inch  using 
Powell's  method,  stress  ratio  method  with  displace¬ 
ment  scaling  and  the  design  algorithm  of  section 
III.  The  final  designs  using  the  three  techniques 
are  presented  in  Table  3.  Since  there  are  no  pre¬ 
viously  published  results  available  for  comparison, 
the  results  of  the  latter  two  methods  are  compared 
with  those  obtained  by  Powell's  method.  The  final 
design  obtained  by  Powell’s  method  is  believed  to  be 
close  to  the  optimum  design,  since  the  method  com¬ 
pared  favorably  with  the  one  by  Schmit  and  Miura  for 
the  2.0  inch  case.  The  design  obtained  by  the  algo¬ 
rithm  of  section  Ill  was  very  close  to  the  Powell's 
method  design.  The  final  design  obtained  by  stress 
ratio  method  with  displacement  scaling  is  consider¬ 
ably  different  with  a  much  higher  weight;  in  fact  it 
has  almost  twice  the  minimum  weight  obtained  for  an 
allowable  displacement  of  2.0  inches  using  the  same 
technique.  It  has  been  observed  in  this  wing  box 
problem  that,  if  a  displacement  limit  of  a  little 
over  2.0  inches  was  imposed,  then  the  optimal  design 
would  be  one  in  which  no  displacement  constraint  was 
active.  Thus  imposing  a  displacement  limit  of  2.0 
inches  or  less  would  result  in  an  active  displace¬ 
ment  constraint  at  the  optimal  design,  but  with  a 
limit  of  over  2.0  inches  no  displacement  constraint 
would  be  active.  Since  the  stress  ratio  method  only 
considers  the  displacement  limit  in  a  crude  way,  by 
simple  scaling,  then  it  might  be  expected  to  have 
difficulty  with  a  problem  where  the  displacement 
constraint  was  critical  at  the  optimal  design,  which 
is  the  case  for  this  example  with  a  displacement 
limit  of  1.0  inch.  On  the  other  hand,  with  a  dis¬ 
placement  limit  of  2.0  inches,  even  though  it  is 
active  at  the  optimal  design,  it  is  not  a  critical 
constraint.  Thus  the  stress  ratio  method  should 
work  well. 

The  final  designs  using  the  three  different 
techniques  for  an  allowable  displacement  of  4.0 
inches  are  presented  in  Table  4.  The  design  obtain¬ 
ed  by  the  algorithm  of  section  III  is  close  to  the 
optimum  design  obtained  by  Powell's  method,  but  is 
slightly  heavier)  however,  most  of  the  design 
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v.u  i.ilil.-;,  .it*-  V'*r y  itiinil.u.  The  upl  imutn  design  ob¬ 
tain'd  t'V  I  Ilf  Sl  M  1 .11  in  llietllilil  Willi  displacement 
scaling  !■;  tw.trly  til* ’  came  .is  Powell':;  method  ex¬ 
cept  the  areas  of  truss  elements  j,  4  are  slightly 
higher.  It  is  interesting  to  note  that  in  this 
method  different  starting  points  made  no  difference 
in  the  number  of  iterations  or  minimum  weight  a- 
chievod.  It  is  also  observed  that  the  stress  ratio 
method  with  displacement  scaling  worked  better  than 
the  design  algorithm  of  section  III  as  it  did  in 
2.0  inch  case.  This  is  due  to  the  fact  that  the 
displacement  constraint  is  not  active  at  the  opti¬ 
mal  design,  thus  it  need  not  be  considered  in  the 
optimization,  which  is  the  case  for  the  stress 
ratio  method. 

In  conclusion,  the  design  algorithm  presented 
in  section  111  gave  good  results  for  all  three 
allowable  displacements  used,  i.e.  2.0  inch,  4.0 
inch,  and  1.0  inch.  Most  of  the  design  variables 
obtained  in  this  method  for  all  throe  cases  were  a 
good  match  with  those  determined  by  Schmit  and 
Miura  and  by  Powell's  method.  Also,  the  constraints 
that  were  active  were  almost  always  identical  to 
those  of  Schmit  and  Miura,  and  Powell's  method.  In 


tin-  design  algorithm  of  section  III,  a  wide  range 
of  starling  points  were  used,  most  of  which  were 
far  from  optimum  design,  nevertheless  the  optimum 
design  was  obtained  in  from  11-13  iterations  for 
all  cases,  while  in  Schmit  and  Miura,  and  Powell's 
method  the  initial  guess  was  chosen  very  close  to 
the  optimum  design.  If  starting  points  had  been 
selected  farther  from  the  optimum  design  it  would 
likely  had  taken  more  analyses  than  indicated  in 
Tables  1-4.  In  the  design  algorithm  of  this  paper, 
a  wide  range  of  relaxation  parameters  n  were  tried. 
It  was  found  that  n  in  the  range  fz„.u  0.05  to  0.15 
gave  the  best  .results  for  all  cases.  The  stress 
ratio  method  worked  well  and  gave  designs  similar 
to  Schmit  and  Miura  and  Powell's  method  when  the 
displacement  constraints  at  the  optimum  design  were 
not  active  (or  when  just  barely  active)  but  it  did 
not  work  when  the  displacement  constraints  were 
critical.  The  design  algorithm  of  section  III  gave 
good  results  regardless  of  whether  the  displacement 
constraint  was  active  or  not,  therefore  use  of  this 
method  is  recommended  over  the  stress  ratio  method 
with  displacement  scaling. 


Table  1  Final  Designs  for  18  Element  Wing  Box,  Allowable  Displacement  2.0  Inches 


Schmit  and  Miura 


Member 

CST  Model  1 

CST  Model  2 

CST  Model  1 

Gallatly 

This 

No. 

SSP 

SSP 

Shear  Web 

&  Berke 

Gallatly* 

Research 

TRUSS 

Aa ( in2 ) 

Aa (in2) 

A^in2) 

AA(in2) 

Ai(in2) 

Ai(in2) 

1 

4.045 

3.151 

2.229 

0.6505 

1.0431 

0.9735 

2 

0.1001 

0.1000 

0.1001 

0.1001 

0.1036 

0.1301 

3 

0.1001 

0 . 1000 

0.1000 

0.2366 

0.3508 

0.2339 

4 

0.1330 

0.2324 

0.3202 

0.2352 

0.3315 

0.3652 

5 

0.1002 

0.1000 

0.1001 

0.1001 

0.1035 

0.1479 

CST 

Tj (in) 

T. (in) 

Tj (in) 

Tj (in) 

Tj (in) 

Tj  (in) 

1,2 

0.08286 

0.08641 

0.1093 

0.1328 

**0.1441 

0.11131 

3,4 

0.05363 

0.05733 

0.05911 

0.0702 

**0.0599 

0.05139 

5 

0.03786 

0.03932 

0.04098 

0.0449 

0.0435 

0.03555 

SSP 

^  (in) 

Tj (in) 

Tj (in) 

Tj (in) 

Tj (in) 

Tj (in) 

1 

0.3636 

0.3851 

0.09345 

0.0876 

0.0876 

0.53075 

2 

0.2236 

0.2152 

0.09437 

0.0889 

0.0895 

0.22626 

3 

0.1310 

0.1361 

0.07687 

0.0808 

0.0664 

0.11120 

4 

0.1156 

0.1004 

0.07293 

0.0768 

0.0553 

0.16274 

5 

0.09166 

0.09113 

0.07570 

0.0815 

0.0537 

0.10408 

6 

0.02000 

0.02000 

0.02001 

0.0200 

0.0219 

0.02070 

7 

0.02000 

0.02000 

0.02001 

0.0200 

0.0215 

0.02070 

8 

0.03096 

0.03090 

0.02804 

0.0337 

0.0256 

0.02994 

Tina] 

Wt.  (lbs) 

402.97 

403.  3r> 

357.32 

387.67 

389.8 

429.38 

Analyst's 

Ntn'di'fl 

9 

ii 

9 

4*** 

193 

13 

•The  original  design  obtained  by  Gallatly  was  scaled  up  so  that  the  triangular  idealization  of 
the  cover  plates  satisfies  stress  constraints. 

••Each  portion  was  modelled  by  a  quadrilateral  clement  in  the  original  work  by  Gallatly. 
•••Subsequent  iterations  gave  heavier  designs. 
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V.  Com  ■  1  union*. 

As  a  result  of  solviin)  sevr. il  example  i.ise s 
it  has  In-eti  shown  that  the  opt  i  mu  1  i  t  y  rr  i  1  erion 
method  developed  by  Khun,  Willmert,  und  Thornton 
can  bo  extended  to  structures  composed  of  different 
element  types.  Starting  from  any  initial  design, 
the  efficient  technique  is  capable  of  locating  a 
near  optimal  design  in  very  few  analyses,  with  each 
iteration  of  the  method  requiring  little  additional 
calculation  beyond  that  required  for  the  analysis. 
That  is,  no  extensive  derivatives  are  required, 
nor  the  determination  of  large  numbers  of  Lagrange 
multipliers. 
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Table  2  Final  Designs  for  18  Element  Wing  Box, 
Allowable  Displacement  2.0  Inches 


Member 

No. 

Powell's 

Method 

Stress  Ratio 
Method  With 
Displacement 
Scaling 

Design 

Algorithm  of 
Section  III 
n=.l5  n=.io 

(in2) 

(in2) 

(in2) 

(in2) 

TRUSS 

1 

3.1768 

3.8082 

0.8629 

0.9735 

2 

0.1017 

0.1211 

0.1037 

0.1301 

3 

0.1016 

0.8209 

0.1785 

0.2339 

4 

0.1016 

0.4745 

0.4291 

0.3652 

5 

0.1021 

0.1056 

0.1277 

0.1479 

(in) 

(in) 

(in) 

(in) 

CST 

1,2 

0.08575 

0.07120 

0.1)434 

0.11131 

3,4 

0.05556 

0.05003 

0.05108 

0.05139 

5 

0.03724 

0.04000 

0.03493 

0.03555 

SSP 

1 

0.44214 

0.46342 

0.51675 

0.53075 

2 

0.21703 

0.24850 

0.22939 

0.22626 

3 

0.12819 

0.16012 

0.11050 

0.11120 

4 

0.11548 

0.13341 

0.16076 

0.16274 

5 

0.09427 

0.09816 

0.10340 

0.10408 

6 

0.02015 

0.02112 

0.02076 

0.02070 

7 

0.02013 

0.02112 

0.02075 

0.02070 

8 

0.02016 

0.03254 

0.02741 

0.02994 

Final 

Weight. 

406.72 

410.50 

429.80 

429.38 

(lbs) 

Analyse? 

Needed 

4605 

9 

13 

13 

Table  3 

Final  Designs  Using  Different  Techniques 
For  Allowable  Displacement  1.0  Inch 

Design 

Stress  Ratio 

Member 

Powell's 

Algorithm  of 

Method  With  Dis- 

No. 

Method 

Section  III 

placement  Scaling 

(in2) 

(in2) 

(in2) 

TRUSS 

i 

2.2897 

1.0110 

7.5117 

2 

0.1100 

0.2993 

0.2)07 

3 

1.1754 

i .  5429 

1.4704 

4 

1 .  i  9  2 1 

0.8385 

0.7881 

5 

0.1102 

0.2402 

0.2107 

(in) 

(in) 

(in) 

CST 

1,2 

0.17404 

0.18478 

0.14092 

3,4 

0.11533 

0.13036 

0.10193 

5 

0.09889 

0.08873 

0.08051 

SSP 

1 

0.20066 

0.20034 

0.97677 

2 

0.16501 

0.15930 

0.49612 

3 

0.25517 

0.21330 

0.33075 

4 

0.18624 

0.18628 

0.26438 

5 

0.22221 

0.17286 

0.19369 

6 

0.02054 

0.02000 

0.04215 

7 

0.02029 

0.02000 

0.04215 

8 

0.07849 

0.06690 

0.065673 

Final 

Weight 

663.96 

669.26 

838.84 

(lbs) 

Analyses 

Needed 

4083 

11 

10 
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Table  4  Final  Designs  Using  Different  Techniques 
For  Allowable  Displacement  4.0  Inches 


Member 

No. 

Powell's 

Method 

Design 

Algorithm  of 
Section  III 

Stress  Ratio 
Method  With  Dis¬ 
placement  Scaling 

(in2) 

(in2) 

(in2) 

TRUSS 

1 

3.5682 

1.0626 

3.6869 

2 

0.1052 

0.1041 

0. >045 

3 

0,1067 

0.1041 

0.6594 

4 

0.1163 

0.1638 

0.3241 

5 

0.1049 

0.1041 

0.1045 

(in) 

(in) 

(in) 

CST 

1,2 

0.08377 

0.11289 

0.06935 

3,4 

0.05234 

0.05423 

0.05143 

5 

0.03816 

0.03601 

0.04024 

SSP 

1 

0.41270 

0.49950 

0.50901 

2 

0.22714 

0.22511 

0.24523 

3 

0.13464 

0.11124 

0.16935 

4 

0.13326 

0.16008 

0.12981 

5 

0.09455 

0.10697 

0.09508 

6 

0.02047 

0.02082 

0.02091 

7 

0.02044 

0.02082 

0.02091 

8 

0.03119 

0.02944 

0.03293 

Final 

Weight 

406.04 

426.72 

417.98 

(lbs) 

Analyses 

4074 

11 

11 

Needed 

